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We study admittance and energy dissipation in an out-of-equlibrium single electron box. The system 
consists of a small metallic island coupled to a massive reservoir via single tunneling junction. The potential 
of electrons in the island is controlled by an additional gate electrode. The energy dissipation is caused by 
an AC gate voltage. The case of a strong Coulomb blockade is considered. We focus on the regime when 
electron coherence can be neglected but quantum fluctuations of charge are strong due to Coulomb interac- 
tion. We obtain the admittance under the specified conditions. It turns out that the energy dissipation rate 
can be expressed via charge relaxation resistance and renormalized gate capacitance even out of equilibrium. 
We suggest the admittance as a tool for a measurement of the bosonic distribution corresponding collective 
excitations in the system. 
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PACS: 73.23.Hk, 73.43.Nq 

It is well-known that the phenomenon of Coulomb 
blockade is an excellent tool for observation of interac- 
tion effects in single electron devices P~H5]. Recently, 
due to the progress in the field of thermoelectricity the 
Coulomb blockade under out-of-equilibrium conditions 
has come into the focus of the theoretical [SrfTT] and ex- 
perimental research [PlrTB] . The simplest mesoscopic 
system displaying Coulomb blockade is a single electron 
box (SEB). The properties of such a system are essen- 
tially affected by electron coherence and interaction. 

The set-up is as follows (see FiglTJ) . Metallic island is 
coupled to an equilibrium electron reservoir (tempera- 
ture T r ) via tunneling junction. The island is also cou- 
pled capacitively to the gate electrode. The potential 
of the island is controlled by the voltage U g of the gate 
electrode. The distribution function of electrons in the 
reservoir is assumed to be equilibrium (Fermi distribu- 
tion) while the one inside the island is arbitrary. 

The physics of the system is governed by the Thou- 
less energy of an island £/rh, its charging energy E c , 
and the mean single-particle level spacing 5. Through- 
out the paper the Thoulcss energy is considered to be 
the largest scale in the problem. This allows us to treat 
the metallic island as a zero dimensional object with 
vanishing internal resistance. The characteristic energy 
(Ed) of electrons inside the island obeys the condition 
5 <C £d "C EcEt^. This implies that characteristic 
energy is high enough to render the system incoherent 
and low enough to keep it strongly correlated due to 
Coulomb interaction [TS]. The dimensionlcss conduc- 
tance of a tunneling junction is small j<1. 

x ' e-mail: yaroslav.rodionov@gmail.com burmi@itp.ac.ru. 



A single electron box does not allow for conductance 
measurements since there is no DC-transport. This way 
an essential dynamic characteristic becomes the set-up 
admittance, which is a current response to an AC-gate 
voltage U g (t) = Uq + Uq cos fit. 

Paper [TB] sparked both theoretical and experimental 
attention to the admittance of such a set-up [T7H22] . As 
it is well-known, the real part of admittance determines 
energy dissipation in an electric circuit. Classically, the 
average energy dissipation rate of a single electron box 
is given as follows 



Wo = n 2 c 2 g R 



R 



— , hn<^gE c , (1) 



where C g denotes the gate capacitance, e - the electron 
charge, and h = 2wh - the Planck constant. Expres- 
sion ([1} presents us with a natural way of extracting 
the resistance of a system from its dissipation rate. 




Fig. [T] The set-up. A SEB is subjected to a constant 
gate voltage Uo- The current through the tunneling 
junction is caused by a weak AC voltage U g (t). 
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Electron coherence and interaction change the clas- 
sical result for the dissipation. The low-temperature 
(T <C 6) coherent regime was pioneered in Ref. [TB]. It 
was shown that the energy dissipation rate Wq can be 
factorized in accordance with its classical appearance 
(PJ but the definition of physical quantities comprising 
it becomes different. Geometrical capacitance C g should 
be substituted by a new observable: mesoscopic capac- 
itance C^. This leads to the establishment of another 
observable: charge relaxation resistance lZ q such that 
R — > lZ q in Eq. (TTJ. Charge relaxation resistance of 
a coherent system differs drastically from its classical 
counterpart. In particular, the charge relaxation resis- 
tance of a single channel junction was predicted to be 
independent of its transmission and equal to h/ (2e 2 ) at 
zero temperature [16j . However Coulomb interaction in 
Ref. [TB] and subsequent works Q2] was accounted for on 
the level of classical equations of motion only. Recently 
the result for quantization of the charge relaxation resis- 
tance in SEB at T <C S has been rigorously derived [23] . 
The admittance in this low temperature regime was in- 
vestigated experimentally by Gabelli et al. [T5] . 

The knowledge of the charge relaxation resistance has 
been extended to a SEB at the transient temperatures 
when thermal fluctuations smear out electron coherence 
but electron-electron interaction is strong. The expres- 
sion for the energy dissipation at this transient temper- 
atures keeps its classical appearance if one substitutes 
the renormalized gate capacitance C g and the charge re- 
laxation resistance lZ q for C g and R respectively [23]. 
Unlike the latter, C g and lZ q have strong temperature 
and gate voltage U g dependance. 

The recent experiment by Persson et al [22] explored 
the energy dissipation rate at these transient tempera- 
tures. The admittance of SEB was measured at fixed 
frequency as a function of pumping amplitude Uq and 
the DC part of gate voltage Uq in a wide range. The 
theoretical analysis of the data in Ref. [22] was carried 
out under assumption of linear response to the AC gate 
voltage: the electrons inside the island were assumed to 
be in the equilibrium with the reservoir. However it has 
not been verified experimentally. It is natural to expect 
that this assumption is violated for the set of data with 
high values of the amplitude Uq. 

Motivated by the experiment [22] we study the ad- 
mittance of a single electron box under the out-of- 
cquilibrium conditions. We consider the linear response 
of a SEB with arbitrary electron distribution function 
in the island to the AC gate voltage. 

A single electron box is described by the Hamiltonian 



where 



H = H a + H c + Ht 



(2) 



(3) 



describes free electrons in the lead and the island, H c 
describes Coulomb interaction of carriers in the island, 
and H t describes the tunneling. Here operators a\ {d} a ) 
create a carrier in the lead (island) . Then the tunneling 
Hamiltonian is 



H t =X + X\ X = J2tk a 



ald c 



(4) 



The charging Hamiltonian of electrons in the box is 
taken in the capacitive form: 

H c = E c (h d - q) 2 . (5) 

Here E c = e 2 /(2C) denotes the charging energy, q = 
CgUg/e the gate charge, and hd is an operator of a par- 
ticle number in the island fid = yj Q d^ a d a . To char- 
acterize the tunneling it is convenient to introduce the 
Hermitean matrix: 



Qaa' 



(2.) 2 [6(eW)6{e$j\ 



1/2 



Y. i ik^k)t ka .. (6) 



The energies e k , £« are accounted from the Fermi level, 
and the delta-functions should be smoothed on the scale 
SE, such that 6 <C SE <C T r ,ed- The classical di- 
mensionless conductance (in units e 2 /h) of the junc- 
tion between a reservoir and the island can be ex- 
pressed as follows g = J2 a 9aa- Therefore, each non- 
zero eigenvalue of g corresponds to the transmittance of 
some 'transport' channel between a reservoir and the is- 
land [25]. The effective dimensionless conductance (g c h) 
of a 'transport' channel and their effective number (A c h) 
are given by 



9cta'9ct'oi 

aa' 



N ch = 



yi 9aa 



y 9aa'9a'a 



(7) 



The dimensionless conductance becomes g = g c hN c h- In 
what follows we will always assume 



g ch < 1, 7V ch > 1, S«l. 



(8) 



Throughout the paper we keep the units such that 
h = e = ks = f except for the final results. 

In the presence of time dependent gate voltage the 
gate charge q in Eq. ([5]) is changed as q C g U g (t)/e. 
The gate voltage is coupled to the operator of particle 
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Here the scaling parameter Z is defined as 



Fig. [5] Charging energy E c h = E c (n — q) 2 as a function 
of gate charge q. 

number inside the island only. Therefore the admit- 
tance of the system (the response of the charge in the 
island to AC part of the gate voltage) is determined 
by autocorrelation function of fluctuating particle num- 
ber: i8(t)([hd(t),hd(0)]}, where 6(t) is Heaviside step- 
function. Due to the presence of strong Coulomb in- 
teraction the behavior of the autocorrelation function is 
non-trivial. It corresponds to collective bosonic modes 
similar to the case of Fermi liquid where the density- 
density correlator is governed by the electron-hole ex- 
citations [26H28] . The latter determines the behavior 
of the autocorrelation function in the absence of the 
Coulomb interaction. In an out-of-cquilibrium regime 
we generally expect the collective mode distribution to 
be different from the distribution of the electron-hole 
excitations. As shown in Ref. jTTj , the collective mode 
distribution coincides with the one for the electron-hole 
excitations even out of equilibrium: 



B tt (r) = 



f[l-F e d (T)F^( T )}de 
J[Ff(r)-F^(r)]de- 



(9) 



Here function F d ' r (T) is given in terms of the Wigncr 
transform /^ ,r (r) of the electron distribution function 
f d ' r {t,t') inside the island/reservoir: F e d,r (r) = 1 — 
2/g' r (r), where a slow time t = (t + t')/2. In the equi- 
librium F^ r = tanh(e/2T r ) and B u = coth(cj/2T r ). 

Results. - We focus on the most interesting case of 
Coulomb peak: the vicinity of a degeneracy point q = 
fc+1/2 where k is an integer. In this parametric regime 
the transport is dominated by the two closest charging 
states [25] (see Fig. [5]) which in the case of g = are 
separated by the Coulomb gap A = 2E c (k + 1/2 — q). 
Due to the presence of the tunneling (finite g) all the ob- 
scrvables, e.g., A, undergo strong renormalization near 
the Coulomb peak. 

For not very high frequencies n <C max{|A|, T r , Ed} 
we obtained the following expression for admittance of 
the SEB 



C a Z^g AO A B_ 



in 



C 4tt 



-in 



gAB_ 



(10) 



Z(X) 



1 



9 
2tt 2 



A 



-1/2 



A 



Bu. 
2lo 



-did, 



(11) 



and g, A are rcnormalizcd tunneling conductance and 
Coulomb gap respectively: 



g = gZ 2 (\), A = AZ 2 (\). 



(12) 



The integral in Eq. (jlip runs over frequencies E c 3> 
\u)\ ^> loq = max{T r , Sd, |A|}. The energy scale ujq de- 
termines the natural scale at which the RG procedure 
has to be stopped [11] . Within logarithmic accuracy we 
find A = IiiE c /loq. 

We stress that our result (TIT)]) is valid for an arbitrary 
electron distribution. To make predictions more con- 
crete we consider the case of quasi-cquilibrium F d = 
tanhe/2Trf, Td > T r as an example. This regime is typ- 
ical for a SEB with the metallic island. It is achieved 
when the energy relaxation rate due to electron-electron 
interaction in the island l/r ee » gS (see e.g., [IT])- 

The real part of admittance (flU)) at fixed n as a func- 
tion of q is shown in Fig. [3] for the out-of-equilibrium 
regime with Td > T r . At fixed C g , C and g the height of 
the maximum is controlcd by the effective temperature 
of electron-hole excitations T e h = lim / ^^ (A/2)i3^ [3D]. 
As it was shown, T r ^ T h ^ Td and T c h ~ Td In 2 
for Td ^> T r [IT] . Therefore, out-of-equilibrium ad- 
mittance is confined within the boundaries Re Gn,T d < 
Re£n < RcQn,T r , where Qn,T d {Qn,T r ) are equilibrium 
admittances at temperatures Td{T r ). 

The dissipative part of the admittance in a SEB has 
been addressed experimentally via radio-frequency re- 
flectometry measurements. The device was exposed to 
a continuous rf-signal |22j . In the experiment the tun- 
neling conductance was estimated to be equal g = 0.5 




fc + l/2 fc + i 1 

Fig. [3j The real part of admittance of the SEB at fixed 
Q as a function of q. We use g = 0.5, Q = 0.02E C and 
Cg/C = 0.24. See text. 
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- C g /C = 0.24, T T = T d = 0.08£ c 

1 Cg/C = 0.32, % = T d = 1.110T/;, 

C s /C = 0.31, T r = 0.054 E c , T d = 0.08E, 
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Fig- [U The dissipative part of admittance of the SEB at 
fixed H as a function of q. Three curves corresponding 
to three different formulae are presented. Dashed line 
corresponds to Eq. (|10|l with Z = 1, A = A, g = g and 
= cothA/2T r . Dotted line is plotted according to 
Eq. (JTDJl with B A = coth A/2T r . Solid line corresponds 
to Eq. (|10fl with non-equilibrium given by Eq. (J9]). 
We use g = 0.5 and SI = 0.02E C . See text. 

such that the SEB was in the strong Coulomb blockade 
regime. We plot the real part of the admittance (|10D at 
fixed SI as a function of q in Fig. 01 There, for a sake of 
comparison, we present ReQn computed i) in the equi- 
librium without taking into account the renormalization 
effects, i.e., with Z = 1 and B^ = cothA/2T r (dashed 
line); ii) in the equilibrium and with the renormaliza- 
tion effects, i.e., with B^ = cothA/2T r (dotted line); 
iii) out of equilibrium with > T r and B^ determined 
by Eq. © (solid line). In all three cases, we use the 
same values of g, E c and SI corresponding to the experi- 
ment [22] • As one can see from Fig. 01 a slight variation 
of ratios C g /C and T r ^/E c allows us to make curves 
for cases i), ii) and iii) indistinguishable. In Ref. [22] 
it is assumed that the electrons inside the island are 
in the equilibrium with the reservoir and the renormal- 
ization effects are not important (case i) above). Val- 
ues of Cg/C and T r /E c are used as fitting parameters. 
The curves presented in Fig. U] however demonstrate a 
more subtle picture. As one can see, the successful fit- 
ting of the experimental data by a 'theoretical' curve 
gives yet no confidence that these assumptions are sat- 
isfied. Therefore, more careful analysis of the experi- 
mental data of Ref. (22] is needed. 

The electron-hole distribution B^ enters admittance 
in a twofold way. The analytical structure of admit- 
tance as a function of external frequency SI is entirely 
determined by B u at uj = —A. The scaling parameter 
Z arising from the renormalization however contains in- 
formation on B^ in wide domain ujq < \uj\ < E c . Ad- 
mittance (|10p can serve as a tool for direct experimental 



measurement of B u . As such can be the measurement of 
a real part of admittance Re Qq at two different driving 
frequencies. Other possibility would be the simultane- 
ous measurement of the real and imaginary parts of Qq 
at a given frequency |19) . Then one can read out AB_^ 
in the entire span of A by tuning the DC gate voltage 
Uq. Measurements of frequency dependence of Re Qq at 
the Coulomb peak (A = 0) provide an access to the ef- 
fective bosonic temperature T h- Thus the admittance 
of a SEB under AC gate voltage can be used as the ther- 
mometer for the electron-hole excitations similar to the 
Coulomb blockade thermometer based on measurements 
of the differential DC conductance in a single electron 
transistor (SET) [T2] . 

The real part of admittance determines energy dissi- 
pation rate: Wo = (C ff /2C)Rc Gn\Un\ 2 - At quasi-static 
regime — > we find that even out of equilibrium the 
energy dissipation rate factorizes into the product of 
well-defined physical observables in full analogy with 
classical expression ([TJ): 

SI 2 , l9 h da' , . 

Wn = ^C 2 g K q \Un\\ K 9 =jf, C 9^g^- ( 13 ) 

The charge relaxation resistance lZ q and the renormal- 
ized gate capacitance C g are related to physical observ- 
ables formally defined as [31] 

g ' = 47rlm v q' = Q + Re ^ (14) 

Here Q = (fid) is the average charge in the island, the 
correlation function K R (t) = i6(t)([X(t),X^(0)]) and 
the limit uj — > is assumed. The physics behind quan- 
tities (fT4"]) can be understood if we turn from a SEB to 
a SET. In the absence of source-drain voltage a SET is 
physically equivalent to the SEB. The quantity g' then 
coincides with the SET conductance |3"2"1[3"3] . The quan- 
tity q' is specific to Coulomb blockade physics and can 
be addressed as the quasi-particle charge [3T| . 

With the help of definitions (fT4"]) we obtained the 
following results in the out-of-cquilibrium regime (for 
.9 « 1) 

g< = - 1 -gAd K \nB_ A , q> = k + \ + (15) 

Equations (|T5|) generalize the results for g 1 [321133] and 
q' |31j derived under the equilibrium conditions. 

Derivation. - Below we describe the main steps of the 
derivation. Further details will be given in [31] . Follow- 
ing Ref. [2S] , we write the Hamiltonian ([2]) in the trun- 
cated Hilbert space of electrons on the island account- 
ing for two charging states: with Q = k and Q = k + 1 
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only (see Fig. [5]). The projected Hamiltonian then takes 
a form of 2 x 2 matrix acting in the space of these two 
charging states: 

H + H f + AS Z + A 2 /AE C 



where 



Ht = tk a a\d a S + H.c. 



(16) 
(17) 



and S z , S^ 1 = S x ±iS v are ordinary (iso)spin 1/2 opera- 
tors. Admittance is proportional to dynamical (iso)spin 
susceptibility Hf (t) = i8(t)([S z {t), S z (0)}) [SJ: 



0n = -ifiC ff nf 



To deal with spin operators out of equilibrium the gen- 
eralization of Abrikosov's pseudo-fermions (PF) tp^, ip a 
is used [35j|36]. Integrating out electrons in the limit 
Ach S> 1 5 we arrive at the following effective action [TT] 

dti>(id t - ^ + 77)^ + § / $(t)iia-il)(t) 

(19) 



2 7" 8 „ 
xIL i j(t,l/)i>(t')'y j (T + i}(t')dtdt'. 



Here the pseudo-fcrmion fields -0, ^ are understood as 
vectors in the tensor product of isospin and Keldysh 
spaces. We inserted the factor exp(r)'ip'ip) with 77 — > —00 
into the density matrix in order to fulfill the constraint 
ip(t)ip(t) = 1. The matrices er z , a± = (a x ±ia y )/2, and 
71 = t x , 72 = To are the Pauli matrices in (iso)spin and 
Keldysh spaces respectively. stands for the matrix: 



n = 



n 



R.A.K 




— iuj(t—t' ) 



de 

27T' 



(20) 

(21) 
(22) 
(23) 



n^(r) = T i | [^(r)-Fj_(r)] 

H^(r)=2iJ(l-F e d (r)F^(r))^ 
The PF dynamical spin susceptibility is given as [34j : 

-r^( ff + w , e , w )[G^ +( 

de 



+r 



if Ait 



(e + uj, e, oj)G : 



(24) 



where the renormalizcd Green's function [24j 
Z(X) 



£ - & ±igT a (e)' 



& = -r) + trA/2, 



r CT ( £ ) = — ( e -f_ ff )[j;- CT 

87T S-» 
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e-£-« 



(25) 



The pseudo-fermion distribution J 7 ^ is not known a pri- 
ori. It is to be determined sclf-consistently from corre- 
sponding kinetic equation. It obeys [11) : 



(26) 



As was shown in [21] all terms of G R G R and G A G A 
type are controlled by renormalization scheme and can 
be discarded. Then, Eq. ([2~4"[) becomes simplified: 



(is) n« p/ H=4-x;%^ 



s 



ujT 



RAR 



w + 2igY a 



(27) 

where r CT = r cr (^ cr ). The vertex function T RAR solves 
the following Dyson equation 



AH( e + w>e>w) = 1 + ^/^ ff 



4 7 2tt 

xlmnf : (B a -CT)r^ fl (e + u; + a;,e + a;,w). (28) 



(29) 



By using Eqs p6 jl -([27 [) and the solution of Eq. ([28 

r RAR {^ + ij,^,u}) 1 w + 2^(r_ CT - r CT ) 



w + 2igr c 



w cj + 2ig(T- a + r CT ) ' 



we obtain expression (|10[) for the admittance. 

The computation of g' and is entangled with the 
computation of K R (see Eq. (fT^jl l. Using the definition 
of K R {t) in terms of the operators A(i), one can obtain 
the following expression [34j : 



'8tt 



"2T 



-ReP^'i?^'-^ 



ilmV u ,{B u ,-B u ,- u ) 
(F e \ u ,- U F £ r )de. (30) 



Here we introduce the transverse spin susceptibility 
V R (t) = i6(t){[S-(t),S + (0)]). Following Eq. (JUJ one 
straight forwardly establishes: 



9 =9 



J ^-ImD R ud„B u , (31) 



q' = Q 



du) 
An I 2n 



(32) 



The average charge in the island is given in terms of the 
average isospin as: Q = k + 1/2 — (S z ). Using the result 
for the transverse spin susceptibility [llj : 



V. 



1 



Z\X) 



B_ A u> + A + i0+ ' 
we obtain results ([T5|) for g' and q'. 



(33) 
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In summary, the paper addresses the admittance and 
energy dissipation in an out-of-equlibrium single elec- 
tron box under strong Coulomb blockade (g <C 1)- We 
deal with the regime when electron coherence can be ne- 
glected but quantum fluctuations of charge are strong 
due to Coulomb interaction. We derived the expression 
for the admittance at frequencies ft <§; max{T r , Ed, | A|}. 
We found that the energy dissipation rate retains its 
universal appearance in the quasi-stationary limit even 
out of equilibrium. It is achieved in terms of specially 
chosen physical observables: the charge relaxation resis- 
tance and the renormalized gate capacitance. We pro- 
pose the admittance as a tool for a measurement of the 
effective bosonic distribution corresponding to electron- 
hole excitations in the system. 
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92474-MHKC and RAS Programs "Quantum Physics 
of Condensed Matter" and "Fundamentals of nanotech- 
nology and nanomaterials" . 
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